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Abstract 

The operator square root of the Laplacian (— A) 1 / 2 can be obtained from the harmonic 
extension problem to the upper half space as the operator that maps the Dirichlet boundary 
condition to the Neumann condition. In this paper we obtain similar characterizations for 
general fractional powers of the Laplacian and other integro-differential operators. From those 
characterizations we derive some properties of these integro-differential equations from purely 
local arguments in the extension problems. 



1 Introduction 

Let us suppose we have a smooth bounded function / : E™ —* K and we solve the extension 
problem 

u(x,0)=f(x) foneK" (1.1) 

Au(x, y)=0 for x G K" and y > (1.2) 

to obtain a smooth bounded function u. It is well known that — u y (x, 0) = (— A) 1 / 2 /(x), and 
therefore we can realize (—A) 1 / 2 as the operator T : f i— » — u y (x, 0) in the above extension 
problem. 

This is easy to show by applying T twice. When we place — u y (x,0) instead of / as the 
Dirichlet condition in (|1. 1|) . we obtain —u y (x,y) instead of u as the solution of (|l.ip - (|1.2p . Then 
T(T{f))(x) = T(-u y (x,0)){x) = u yy (x,0) = -A x f{x). To show T = (-A) 1 / 2 it is only left to 
check that T is indeed a positive operator, which follows by a simple integration by parts argument. 

In this paper we will generalize this situation to a similar extension problem for each fractional 
power of the Laplacian. We will construct any fractional Laplacian from an extension problem 
to the upper half space for a specific elliptic partial differential equation. This allows us (in 
forthcoming work) to treat non linear variational problems, involving fractional Laplacians, with 
standard local perturbation methods from the calculus of variations. Examples of such methods 
are the Almgren Frequency formula, and the Boundary Harnack inequality presented below. The 
partial differential equation that we obtain will turn out to be degenerate for any power of the 
Laplacian other than (—A) 1 / 2 , however they belong to a more general class of equations that 
shares many of the essential properties of uniformly elliptic equations (as in [7], [BJ, [5] for the 
divergence case, or [5] for the non-divergence case). 

The fractional Laplacian of a function / : K™ — > M. is expressed by the formula 

(-A)'/(s) = C n , s [ f (x) -J12 d£ (1.3) 
where the parameter s is a real number between and 1, and C n s is some normalization constant. 



It can also be denned as a pseudo-differential operator 

(^AP/(0 = ier/(£) 

The fractional Laplacian can be defined in a distributional sense for functions that are not 
differentiable as long as / is not too singular at the origin or, in terms of the x variable, as long as 

ax < +00 



+2s 



(i + N)» 

We will relate the fractional Laplacian with solutions of the following extension problem. For 
a function / : 1™ — >■ R, we consider the extension u : R™ x [0, 00) — >■ R that satisfies the equation 



u(x,0) = f(x) (1.4) 



a 



Aj;KH Uy + Uyy = (l'^) 

The equation (|1.5p can also be written as 

div(j/ a Vu) = (1.6) 
Which is clearly the Euler-Lagrange equation for the functional 

J(u) = J \Vu\ 2 y a dX (1.7) 

y>0 

We will show that 

an«,. „ 1 u(x. y) — u(x, 0) 
C(—A) f = lim -y a u v = lim ^ ,y \ 

for s = ijp and some constant C depending on n and s. Which reduces to the regular normal 
derivative in the case a = (as in (|l.l[) - Q1.2[l ). 

/ \ i — a 

If we make the change of variables z = I ) in (|1.5p , we obtain a nondivergence form 
equation 

A x u + z a u zz = (1.8) 

for a = Moreover, y a M y = u z . Thus, we claim that the following equality holds up to a 

multiplicative constant 

(-A) s f(x) = - lim . y a u y (x,y) = -u z (x,0) 

8/->0+ 

It seems convenient to keep this notation. The variable x is always in R™. The variables y 
and z are nonnegative real numbers that satisfy the relation z = (i~z<jl • ^ nc f unc tion u is 
the extension of / to the upper half space and takes values in R n x [0, oo). We will use u and its 
version after the change of variables y >— > z indistinctly, and we will call the variable either y or z 
to point out the difference. Whenever we refer to a point in R™ +1 we will use capital letters (like 
X). 



2 Properties of the PDEs 

In this section we will study the basic properties of the equations (II. 5|) and (|1.8p . We will develop 
explicit Poisson formulas among other properties. 
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2.1 Harmonic functions in n + 1 + a dimensions 

The equation (|1.5p has a curious intuition behind it that will help us obtain several properties of 
its solutions. 

For a nonnegative integer a, suppose u{x, y) : 1" x R 1+a — > R is radially symmetric in the y 
variable, meaning that if \y\ = \y'\ = then u(x, y) = u{x, y'). We can think of it as a function of 
x and r, and in these variables write an expresion for its Laplacian: 

a 

An = A x u H — u r + u rr 
r 

We thus obtain an identical expression of equation p.5[) (with r instead of y, but that's just 
naming). As far as the expresion is concerned, there is no need to consider only integer values 
of a. We can then think of the equation (|1.5[) as the harmonic extension problem of / in 1 + a 
dimensions more. Although it is impossible to think of a meaning for jj n + 1 + a when a is not an 
integer, the solutions of (jl.5[) will satisfy many properties common to harmonic functions. 



2.2 Fundamental solution at the origin 

To obtain the fundamental solution for (|1.5[) at the origin, we just have to consider the fundamental 
solution of the Laplacian in n + 1 + a dimensions. If n — 1 + a > 1, it reads 

T(X) = C n+1+a — (2.1) 

where the constant C n +i+ a is given by the formula C k = TT k ^T(k/2 - l)/4. 

It can be verified as a direct computation that T is a solution of (|1.5p when y ^ and indeed 
limj,_^ + y au y — ~CSq for some constant C as we will see later. Notice that r(x,0) = . y£=rrz is 

the fundamental solution of the fractional Laplacian (— A)~r~ for some appropriate constant C 
depending on n and a (Recall X = (x,y)). 

Since the equation (JT78J) can be derived from (|1.5|) by just a change of variables, we can also 
change variables to obtain a corresponding fundamental solution 

f(x,z) = C n+1+a ^t+t (2.2) 

(N 2 + (l-a) 2 |z| 2 ^) 2 

that solves (jl.8p when z ^ and u z (x, z) — > — 5q as z — > 0. 



2.3 Conjugate equation 

We have seen that if u is a solution of (jl.5p . then the function w(x, y) :— y a u v {x, y) seems to carry 
interesting information. Indeed, this function satisfies the conjugate equation 

a 

A X W - -Wy + Wyy = 

that is nothing but (|1.5|) with —a instead of a. This property of w can be checked by a direct 
computation: 

a ( \ a? a a(a — 1) a 

A X W - -Wy + Wyy = J/ [<JyA X U - ~Uy - -Uyy H ~ Uy + 2~Uyy + Uyyy 

a fa A a a 

V I (JyA x ii ^yy ' ^yyy 

= y a 8y ^A X U + \ly + Uy 
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The intuition behind the above fact is simple to explain when n = 1 (i£l and y G [0, +00)). 
The function w turns out to be the stream function related to u x in the following sense: if we set 
v = u Xl then clearly v is a solution of (|1.5|) . since the equation is invariant under translations in 
x. Thus we have div(y a \/v) = 0. The vector field (y a v y , —y a v x ) is then irrotational and there is 
a function w such that Vu> = (—y a v y ,y a v x ). This function w is the stream function of v, and it 
satisfies the equation div(y _a Vw) = 0. Now we check that w = y a u y (x, y), since 

= (y a v v , -y a v x ) 



a 




~Uy - 




\y 





2.4 Poisson formula 

We want to develop a Poisson formula P to explicitly solve (|1.4[) - p.5[) . 

u(X) = J P(x-£,y)f(0 d£ 

The Poisson kernel P must be a solution to (|1.5[) where y > and lim^o P(%, y) — <^o- from 
the previous sections, we see that the correct choice would be to take P(x, y) = —y~ a d y T_ a (x, y). 
Therefore ^ 

P{x,y) = C n>a V - „ +1 _ a (2.3) 

(N 2 + |y| 2 ) 2 

The fact that P is a solution to (|1.5|) where y > is a consequence of the fact that r^ Q is 
to the conjugate equation. Of course it can also be checked by a direct computation. Moreover, 
P(x,y) — y~ n P(x/y, 1), so P must converge to a multiple of the Dirac delta as y — > 0, so P is 
indeed the Poisson kernel. 

The corresponding Poisson kernel for the equation (JTTHJ) can be obtained either by a change of 
variables from (|2.3p or by computing —P z : 

P(x, z) = C n , a (2.4) 

(M 2 + (l-a) 2 M 2/(1 - Q) ) 2 



3 Relation with fractional Laplacian 

In this section we will see how the equations (|1.5[) or (|1.8p relate to the operator (— A) s . Namely, 
we will show that up to a constant factor 

\im y a u y (x,y) = u z (x,0) = -(-A) s /W = / ff° ~.[if s d£ (3.1) 

J \t — x\ 
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3.1 Proof using the Poisson formula 

We can compute u z (x,0) using the Poisson formula. 

, u(x, z) — u(x, 0) ,„ „. 

uJx, 0) = lim v ' ' (3.2) 

2^0 z 

= lim - f P(x - £, z)(f(0 - f{x)) d£ (3.3) 

z^O z J 

r 

= lim / 5n=i-(/(0 ~ /(*)) d£ (3.4) 

i (irr-^ + a-a^lzl 2 ^ 



^PV / d f (3.5) 



= -CC-A)-^/^) (3.6) 

where the limit in (|3.4[1 exists and equals ()3.5() as long as / is regular enough. 

Changing variables, this implies that lim y ^o ~y au y( x iy) a ls° converges to a multiple of the 
fractional Laplacian. The fact that y a u y (x 1 y) has a limit as y — > 0, immediately implies that 
lirxiy-kO "^'^1-"^'°^ has the same limit. 

Remark 3.1. The proof was done by computing u z . If instead we had chosen to compute y a u y 
or u ( x 'V)-_u(. x fi) com plexity of the computation would be comparable. 



3.2 Proof using Fourier transform 

Alternatively, we can also prove (|3.ip taking Fourier transform in x. 

One way to do this is proving that the corresponding energy functionals coincide. Namely 



|Vu| 2 y a dX 



The equation (|1.5p becomes 

- I£l 2 u(t y) + y) + Uyy(£, y) = o 

We thus obtain an ordinary differential equation for each value of £. 
Suppose that <f> : [0, oo) — > M is the minimizcr of the functional 



(3.7) 



J(<f>) := y (|^r + l^l 2 ) 2/ Q d?/ 



for (f)(0) = 1, then solves the following equation: 



-<j)(y) + ~<f>y(y) + (j>yy(y) = 

0(0) = 1 
lim 4>(y) = 

y— >oo 
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By a simple scaling we can see that 

Ht,y) = f(0<t>(\t\y) 

and thus, the energy of u becomes 



J \Vu\ 2 y a dX = J Jm 2 \u\ 2 + \u y \ 2 )y a dydC 

/(Ol 2 |e| 2 fl0(ie|y)| 2 + |^'(|e|y)r)y a dy d£ 



\m\ 2 i£i 1_0 / (\m\ 2 + wm 2 )y a ^ 



|/(0| 2 |^| 1_O ^) d£ 

Thus we conclude (I3.7[) . The corresponding Euler Lagrange equations for each energy must 
then coincide up to a constant factor. Therefore we obtain 

- \imy a u y {x,y) = C(-A) 1 ^ f(x) 

y—>0 

where the constant C is given by J \4>). 

We could also make another proof from (jl.8[) . After taking Fourier transform, the equation 
(|1.8p becomes 

-\H\ 2 u(ti,z) + z a u zz (tz) = 

We thus obtain an ordinary differential equation for each value of £. 
Suppose that <f> : [0, oo) — > M solves the following equation 

-z a <j>"(z) + <f>(z) = 

0(0) = 1 (3.8) 
lim 4>(z) = 

z — >oo 

Then by a simple scaling we can see that 
And then, 

fi,(£,0)=fi(£,0)|£|^'(0) 

= c a |^| 1 - a ^,o) = c a |^| 1 - a /(0 

To prove (|3.ip we now only need to show that such function <f> exists and it is differentiable at 
zero. We notice that for e small enough then 

4>{z) := min(l, z~ £ ) 

is a supersolution, whereas for A and B large 

cp{z) := e -^ 1/2 +^ 2 



G 



is a supersolution. Thus, by Perron's method, we can find a <p in between that solves 
Moreover this cj> is Lipschitz at z — and < <f> < 1 since it is trapped between <fi and 0. From the 
equation (|3.8p we also see that then < <j) zz < z~ a , thus S C 2 ~ Q . Not only does this complete 
the proof of (|3.ip but it also gives an interesting regularity estimate for the solution u with respect 
to z. In the next section we will explore this issue further. 

4 Reflection extensions 

In order to apply interior Harnack estimates to our PDE (|1.5[) (or (|1.8p ). we must show that if the 
operator (— A) s / = in an open ball, then we can reflect the solution u and make it a solution of 
(|1.5[) (or (|1.8[) ) across y = (z = 0) in a suitable sense. 
We first address the divergence case. 

Lemma 4.1. Suppose that u : R" x [0, oo) — > R is a solution to (|1.5p si/c/i i/ia£ /or |x| < r 

limj\(x,|/)=0 (4.1) 



£/ien i/ie extension to the whole space 

u(x,y) 

is a solution to 



u(x,y) y>0 
u(x,-y) y<0 



div(|w| a Vu) = 

in the week sense in the (n + 1) dimensional ball of radius R ({(x,y) : \x\ 2 + \y\ 2 < R 2 })- 
Proof. Let h G C§°(Br) be a test function. We want to show that 



J Vu-Vh \y\ a dX = (4.2) 



B r 

Let e > 0, we separate a strip of width e around y = in the domain of the integral above 



J X7u-X7h\y\ a dX = J Vu ■ Vh \y\ a dX + J Vu ■ Vh \y\ a dX 

Br B R \{\y\<s} B R n{\y\<e} 

J div{\y\ a hVu) dX+ J Vu-Vh\y\ a dX 

Br\{\v\<s} B a n{\y\<e} 

J hu y (x,e) e a dx+ J Vu ■ Vh \y\ a dX 



B R n{\y\=e} B a n{\y\<e} 

When we let e — > 0, the second term of the right hand side above clearly goes to zero because 
\y\ a \Vu\ 2 is locally integrable, and the first term converges to zero if e a u y (x, e) — > as e — » 0. 

Therefore, if e a u y (x,e) ~> as e — > 0, then div(|u| a Vu) = in the whole ball Br (accross 
y = 0). □ 



We must clarify in what sense we take the limit in (|4.ip . In this particular case, since the limit 
vanishes, we could prove that actually u is C°° in x and the limit holds in the uniform sense. In 
general it would be convenient to have a weak definition of (|4.1jl . The equation (|4.2p becomes the 
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definition of (|4.ip in the weak sense. In any case the limit (|4.1[) holds in the sense that, if we take 
any smooth test function tp G Cq°(B*), 

lim / y a Uy(x,y)ip(x) dx = 
v->oj 

B* 

For a general class of functions g, the problem 

div(|j/| a Vm) = in the weak sense in Br 

u = g on OBr 

has a unique solution in the weighted Sobolev space H 1 ' 2 (Br, \y\ a ). This kind of problems and 
general properties of elliptic equations with A 2 weights are studied in [7] . 
We now address the nondivergence case. 

Lemma 4.2. Given a continuous function g on dBji, such that g(x, z) = 17(2;, —z), there exists a 
unique function u G C(Br) such that: 

i. u solves u xx + \z\ a u zz = in Br U {z ^ 0} in the classical sense. 

ii. u G C X (B R ) 

in. u z (x, 0) = 

Moreover, for this solution the Harnack inequality result of J5J/ applies. 

Proof. We point out that any viscosity solution of (|1.8p is C 2 away from z = 0, so it would be a 
solution in the classical sense in Br U {z ^ 0}. 

Let us prove uniqueness first. Let us suppose there were two solutions u and v satisfying i., ii., 
and iii.. For an arbitrary e > 0, consider w = u — v + e\ z\. Then w < eR on 8Br. Let us suppose 
that w has an interior maximum at a point x G Br. This point cannot be in Br fl{z^ 0}, since 
there w solves the non-degenerate elliptic equation (|1.8p and thus it does not have an interior 
maximum. Therefore, if there is any interior maximum, it has to be on {z = 0} n Br. But clearly 
there cannot be a maximum since d^w > d~w. Therefore w < eR in the whole ball Br. Since 
e is arbitrary, we conclude that u < v in Br. Similarly we can obtain that v < u in Br, so they 
must coincide. Notice that we only used (i.) and (ii.) for uniqueness. 

Now, let us prove existence. The subtle point here is to show that there is a C 1 solution. What 
we do is to prove a uniform C 1 '" estimate for solutions to the problem 

u £ = g on dB R (4.3) 

A x u e + (\z\ + e) a u e zz = in dB R (4.4) 

For any e > 0, by the Schauder theory, this problem has classical solutions. If we have a C 1 ^ 
estimate uniform in e, we take limit as e — > and obtain the desired solution. 

The solutions u £ are uniformly bounded in L°° due to the maximum principle for uniformilly 
elliptic equations. The equation (|4.4[) has constant coefficients with respect to x. This means that 
we can differentiate the equation in a standard way to obtain uniform estimates for derivatives of 
any order with respect to x (In terms of fractional Laplacian, this corresponds to the fact that 
functions such that (— A) s u = are C°°). Therefore we have that A x u e is bounded independently 
of e in any smaller ball B^ 1 _g/ 2 )R- 

Since u £ G C 2 (Br) and u is symmetric respect to the hyperplane z = 0, then u e z (x, 0) = 0. 
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From the equation (I4.4|) : 



< 



lA^I 
(\z\+sY 

c 



Recall that a — and since a € (—1, 1), then a < 1. We can then integrate ii| z for any x, z 
such that \x\ < 1 — 25 and < z < 1 — 5. 



J u zzi x ' s ) ds 



z 



This shows u e y is C" 7 in B(^_^ R for 77 = min(l, 1 — a), independently of e for any 5. 

So, we can take e — > and extract a subsequence that converges to the desired solution u. 

Moreover, for any e, u £ is smooth and satisfies the equation (|4.4[) that has smooth coefficients 
and for which the Harnack inequality from [5] applies. Therefore, the same estimate passes to the 
limit as e — > and the solution u of the original problem satisfies Harnack inequality. □ 

Remark 4.3. We have to be careful if we want to study the viscosity solutions to (|1.8|) . The naive 
definition using C 2 test functions would not suffice for uniqueness. If, for example, a > and we 
consider the function u(x, z) = \z\, then for any C 2 function <j> touching u from below at a point 
(x, 0), A x <p(x,0) + a 4>zz(x,0) — A x cj) < and thus u(x,z) = \z\ would be a (non-differentiable) 
viscosity solution. However if we test against less regular functions, like <f>(x,z) = z 2 ~ a , then u 
would fail to be a supersolution. In the notation of [4], this corresponds to the distinction between 
C-viscosity solutions and L p -viscosity solutions. 



1*4(^2)1 



5 Harnack and boundary Harnack type estimates 

As an application on how to apply the equations (|1.5p or (|1.8p to the study of fractional harmonic 
functions, we prove a Harnack inequality and a boundary Harnack inequality using (local) pde 
methods. 

Harnack inequality is not a new result for the fractional Laplacian. It can actually be proved 
using direct classical potential methods like in ^9 . A boundary Harnack estimate for the fractional 
Laplacian was first proved in [3] using potential methods. 

In this paper we derive the Harnack and boundary Harnack inequality for the fractional Lapla- 
cian from the Harnack inequality for singular elliptic equations, either with A2 weights (see [7], 
[H] and also [TU]) or for certain classes of nondivergence problems (see [5]). 

Theorem 5.1 (Harnack inequality). Let f : R" — ► R be nonnegative such that (—A) s f — in 
B r . Then there is a constant C (depending only on s and dimension) such that 

sup / < C inf / 

B r /2 B r/2 
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Proof using (|1.5|) . We consider the extension u of / that solves (II . 5[) . This function u is going to 
be nonnegative since / is. We reflect it through the {y = 0} hyperplane. Since (—A) s f — in 
B r , by Lemma |4. 11 u is a solution to 

div {\y\ a Vu) = 

in the n + 1 dimensional ball of radius r centered at the origin. 

We can then apply the result of [7] to obtain Harnack inequality for u and thus also for / in 
half of the ball. □ 

Proof using (|1.8j) . We consider the extension u of / that solves (|1.8|> . This function u is going to 
be nonnegative since / is. We reflect it through the {z = 0} hyperplane. Since (— A) s / = in 
B r , then u satisfies the conditions of Lemma [4.21 We can then apply the result of [5] to obtain 
Harnack inequality for u and thus also for / in half of the ball. □ 

Remark 5.2. In the two proofs above we can observe that what is needed for Harnack inequality 
is that the function u is nonnegative in an n + 1 dimensional ball. The condition / > is indeed 
sufficient for that, but it is not strictly necessary. 



Theorem 5.3 (Boundary Harnack). Let f,g : R™ — > R be two nonnegative functions such that 
(— A) s / = (— A) s g = in a domain tt. Suppose that xq £ d£l, f(x) = g(x) = for any x G B\\VL, 
and dfl n B\ is a Lipschitz graph in the direction of x± with Lipschitz constant less than 1. Then 
there is a constant C depending only on dimension such that 

sup 44 <C inf 44 (5-1) 
xennsi 9{x) zefinBi g(x) 
u 2 



for any x, y G tt (~) B 1/2 (x ). 

Proof. We consider the extension ttW of / and u' 2 ' of g that solve (|1.5p . As before, we reflect 
through {y = 0} for k = 1, 2 and obtain a solution across this hyperplane through f2. What we 
want to do is to find a map that straightens up the domain in order to apply the result of [8]. 

First, since fi is a Lipschitz domain, we can find a bilipschitz map ipi : R" — > R" such that 
^1(^0) = and <pi(f2) n -B1/2 = B1/2 H {xi > 0}. We can extend this map to R ra+1 as constant in 
the variable y. Now, the functions — 'Pi 1 are solutions of the equation 

div (\y\ a bijdjU^) =0 in {y 7^ 0} and across {y = 0} on (fii(tt) 

where the matrix 6y is given by ^jj— ■ Since <pi is bilipschitz, then the singular values of ZJt^i 
are bounded below and above and 6^ is uniformilly elliptic. 

Now we take the map tp 2 that maps M™ +1 \{xi < OAy = 0} into the half space R n+1 n{a;i > 0}. 
This map is constant in the variables x\, . . . , x n , and if we write the pair {xi,y) in polar coordinates 
(r, 9), it maps it to (r, 8/2). The singular values ofDtf2 are all equal to 1 but the one in the direction 
of ^ ( i.e. the one in the direction (— y, 0, . . . , 0, Xi)), that is 1/2. Thus the singular values of 

(k) (k) —l 

Difi2 are bounded above and below. The functions u 3 — u 2 o ip 2 satisfy the equation 

div(cy5 J 4' C) ) = in B i/2 n {xi > 0} 

The matrix Cij is given by (cy) = Dip 2 (bij) D(p 2 m(X), where the weighted measure m(X) dX is 
the pull back of the measure \y\ a dX by <p% . Since the A 2 class of weights is invariant by (local) 
bilipschitz transformations, the equation still fits into the theory in [8]. 
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We can apply the boundary Harnack inequality of [8j to the functions u^' and to obtain 

sup f < C mf 



Bin{xi>0} u { p{x) Bin{x x >0} u f\x) 



Recalling that f(x) — o ip 2 o . . . , x n ,0) and g(x) = U3 o ip 2 o ipi(x\, . . . ,x n ,0), we 

obtain that 

sup 44 <C inf 44 (5-2) 
leans, g{x) xennB r g[x) 



for some universal r < 1. By a standard covering argument, (|5.2[) together with Theorem 15.1 
imply (|5~Tj) . □ 



As before, the condition in Theorem 15.31 that requires / and g to be nonnegative in the whole 
K ra is not sharp. In the proof we can see that what we need to apply the result of [8] is that 
the corresponding extensions vS 1 ^ and are nonnegative in the unit ball of n + 1 dimensions. 
The advantage is that this is a local condition, therefore we can obtain the standard corollary of 

(!)... . f 

boundary Harnack saying that is C a in a neighborhood of the origin. In particular ^ is C a , 
which is not strictly a corollary of Theorem 15.31 but of its proof. 



Corollary 5.4. Let f and g be as in Theorem 1 5. 31 then is a C a junction in -B1/2 H f2 for 
some universal a 

Remark 5.5. In [8], the boundary Harnack principle is proven for A2 weights and Lipschitz 
domains. An extension of that result to nontangentiably accesible (NTA) domains (as in the last 
section of [2]) would lead to a boudary Harnack principle for the fractional Laplacian to domains 
only requiring a uniform capacity condition. 



6 Monotonicity Formulas 

Monotonicity formulas are a very powerfull tool in the study of the regularity properties of elliptic 
PDEs. They have been used in a number of problems to exploit the local properties of the equations 
by giving information about the blowup configurations. Since the equation (|1.5p represents a 
harmonic function in n + 1 + a dimensions, it is to be expected that any monotonicity formula 
known for harmonic functions will have its counterpart for solutions to (|1.5| . For example, the 
simplest one would be that if 

A x u H — u y + u yy = for (x, y) E 

lim y a u y (x,y) = Ofor \x\ < 1 

y->0 



then 



1 f 2 

^(R) = j^ n+ i +a I \ ^ u \ J/ a dA is monotone nondecreasing in R 



where stands for the n + 1 dimensional half ball {(x, y) : \x\ 2 + y 2 < 1 A y > 0}. 

The weight y a is the correct one if we think that the function u is radially symmetric in a + 1 
variables and the y is the modulus of those variables as it is explained in section 12.11 

More interestingly, we have Almgren's frequency formula 
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Theorem 6.1. If u is a solution to (| 1 . 5|) in B^ o such that for any x in B^ a , either u(x,0) = 
or liniy^o y au y{x, y) = 0, then 

f B+ \Vu\ 2 y a dX 

<&(-/?) = R — - g * s monotone nondecreasing in R for R < Rq 

JdB+ M y ada 

Moreover, $ is constant if and only if u is homogeneous. 

The proof is essentially the same as for harmonic functions and it is mainly computational. 
The original proof can be found in pQ. Here we only need a few minor modifications for our case. 
First we need the following lemma. 

Lemma 6.2. Let u be s solution of (|1.5p in Br n {y 0}, such that y a u y (x,y) is bounded and 
for every x with \x\ < Rq either u(x, 0) = or lim^o y a u y (x, y) = 0. Then the following identity 
holds for any R < Rq. 

R J (K| 2 -K| 2 ) y a da= J (n + a-l)\Vu\ 2 y a dx (6.1) 

dB a n{y>0} B+ 

where u T stands for the gradient of u tangential to dB r . 

Proof. The only thing we have to notice is that since u solves (|1.5p . then 

div (y A ^fx - y"(X, Vu)Wu^ = y"^±^± |Vu| 2 

So we apply divergence theorem in the set Br n {y > e} to obtain 

R I y a (^-W 2 ) da+ f _ £ ° I V "^' £ )I e + e a ((x,s),\7u(x 1 e))u y (x 1 e)dx 



dB R n{y>e} x ' \x\<R 



n + a — 1 „ ._ l2 , 
y a \Vu\ 2 dx 



The second term comes from the integration of the vector field at the bottom of the half sphere. 
We observe it goes to zero as e — > since there we have that for each x either u(x,0) = or 
lim y ^Q y a u y (x,y) — and a < 1. So we can remove that term and extend the expresion to the 
whole ball Br. 

I Vm ! 2 . |2 \ a i f 71 + a - 1 ir, |2 



R J I — \ u v\ \ y a da= J — — |v«r y a dx 

dB R n{y>0} ^ ' B j 



Expanding |Vu| 2 = \u T \ 2 + \u v \ z we obtain (|6.ip . □ 



Proof of Theorem \6.1\ We will show that log $(-??,) is nondecreasing. By the usual scaling argu- 
ment, it is enough to check it at r = 1. We compute log<i>: 

log$(i?) = logi?+log / |Vu| 2 y a dX - log J \u\ 2 y a d9 



B+ dB + 



12 



and now we compute its derivative at R = 1. Let Si = dB\ f) {y > 0}, then 

L |Vu| 2 y a da L (2uu v + (n + a) \u\ 2 ) y a da 

(log$)'(l) = l + -^ " y > 1 1 ,y (6 . 2) 

J Bt \Vu\ 2 y a dX J Si \u\ 2 y"da 

Now we observe that since div(y a Vu) = 0, we have div y a uVu = y a |Vit| 2 , and then 

J |Vu| 2 y a dX = lim + J uu u y a da + J u y uy a dx 



Si \x\<l 

uu v y a da 



Si 

From Lemma 16.21 we have that 



J \Vu\ 2 y a da = J (K| 2 -K| 2 ) y a da + 2 J \u v \ 2 y a da 

Si Si Si 

= J (?i + a-l)|Vu| 2 y a dX + 2 J \u u \ 2 y a da 



B? 



Putting (fOj) and (|6T4|) together with (|6~2)) we obtain 



(6.3) 



(6.4) 



J R+ (n + a - 1) |V«r y a dX + 2 J \u v \ V a &° 

(log $)'(!) = 1 + , * 

J 5l mw. y a da 

J s (2uu v + (n + a) \u\ 2 ) y a der 

, 2/ s J^| 2 fd. j 2uu v y°da 

= 1 + (n + a — 1) — in + a) H ? ; * 

/ 5l uu » y d(7 \ u \ 2 y»da 

/ J Si K| 2 j/ a dg ^ y a da \ 

V yad(T Isi M 2 y ada ) 

Therefore (log$)'(l) > follows from Cauchy-Schwartz inequality since 
J uu v y a da j < ( J \u\ 2 y a da J A J \u„\ 2 y a da 

\Si / \Si I \Si / 

and the equality is achieved if and only if u v = Xu for some constant A. Therefore, only when u 
is homogeneous of degree A. □ 

7 Other operators 

The results in this paper suggest that many integro-differential operators can be thought in the 
same way. Let us suppose we have smooth uniformilly elliptic coefficients a%j in K n x [0, oo). 
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Consider the following extension problem: 

u(x,0) = f(x) for i e l n (7.1) 

a i:j u i: j (x. y) =0 for x G R" and y > (7.2) 

If we now consider the operator / i— » u y (x, 0), this is going to be some integro-differential operator 
of degree one that we can also obtain from a local type pde. We can push the situation even 
further by considering coefficients ciij that are only measurable, and maybe singular near y = 0. 
It is not reasonable to expect that every integro-differential operator can be realized in this way. 
It looks difficult to characterize the ones that can. But certainly for many it is possible. 

Let us consider the case when the operators are invariant under translation in x. We want to 
study the operator that maps the Dirichlet to the Neumann condition for the equation 

A x u + a(z)u zz = 

we take the Fourier transform in x and proceed as in section [321 to obtain the differential equation 

-aW(z) + |e| 2 <Xz) = 

m = i (7.3) 

lim <f>(z) = 

z — >oo 

The operator T : u(x, 0) i— » —u z (x, 0) is then the pseudodifferential operator whose symbol s(£) is 
given by solving (|7.3p and computing <j) z (0) for each value of £. 

The question is what symbols s(£) can we obtain by the above procedure. We can see that 
s(£) is radially symmetric and monotone increasing in |£|. Indeed, if <f>\ is a solution of (|7.3[) for 
|£| = ro and ipi for |£i| = r\ > tq then 

-a{z)<p'{{z) + \!io\ 2 Mz) = (ICol 2 - \tif)Mz) < 

Thus 0i < 0o by comparison principle. 

As £ — > 0, the solution <f> of (|7.3p converges to <fi = 1, so s(0) = 0. We leave the question if 
every radially symmetric symbol s such that s(0) = and it is monotone increasing with respect 
to |£| can be obtained by a suitable choice of a(z). 
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